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A PIEZOELECTRIC EULER-BERNOULLI BEAM WITH DYNAMIC 
BOUNDARY CONTROL: STABILITY AND DISSIPATIVE FEM 

MAJA MILETIC AND ANTON ARNOLD 


Abstract. We present a mathematical and numerical analysis on a control model for 
the time evolution of a multi-layered piezoelectric cantilever with tip mass and moment 
of inertia, as developed by Kugi and Thull m- This closed-loop control system consists 
of the inhomogeneous Euler-Bernoulli beam equation coupled to an ODE system that 
is designed to track both the position and angle of the tip mass for a given reference 
trajectory. This dynamic controller only employs first order spatial derivatives, in order 
to make the system technically realizable with piezoelectric sensors. From the literature 
it is known that it is asymptotically stable m- But in a refined analysis we first prove 
that this system is not exponentially stable. 

In the second part of this paper, we construct a dissipative finite element method, 
based on piecewise cubic Hermitian shape functions and a Crank-Nicolson time discretiza¬ 
tion. For both the spatial semi-discretization and the full x — t-discretization we prove 
that the numerical method is structure preserving, i.e. it dissipates energy, analogous to 
the continuous case. Finally, we derive error bounds for both cases and illustrate the 
predicted convergence rates in a simulation example. 


1. Model 

The Euler-Bernoulli beam (EBB) equation with tip mass is a well-established model with 
a wide range of applications: for oscillations in telecommunication antennas, or satellites 
with flexible appendages Eng, flexible wings of micro air vehicles [S], and even vibrations 
of tall buildings due to external forces m- The interest of engineers and mathematicians in 
the corresponding control problems started in the 1980s. So various boundary control laws 
have been devised and mathematically analyzed in the literature - with the stabilization 
of the system being a key objective (cf. [ 31 ]). Soon afterwards, also exponentially stable 
controllers were developed which require, however, higher order boundary controls for an 
EBB with both applied tip mass and moment of inertia m- On the other hand, if only a tip 
mass is applied, lower order controls are sufficient for exponential stabilization [12] . In spite 
of this progress, and due to its widespread technological applications, considerable research 
on EBB-control problems is still underway: In the more recent papers jH] [30] exponential 
stability of related control systems was established by verifying the Riesz basis property. 
For the exponential stability of a more general class of boundary control systems (including 
the Timoshenko beam) in the port-Hamiltonian approach we refer to [49j . 

We shall analyze an inhomogeneous multi-layered piezoelectric EBB with applied tip mass 
and moment of inertia, coupled to a dynamic controller that uses only low order boundary 
measurements. This system was introduced by Kugi and Thull in m to independently 
control the tip position and the tip angle of a piezoelectric cantilever along prescribed tra¬ 
jectories. This beam is composed of piezoelectric layers and the electrode shape of the layers 
was used as an additional degree of freedom in the controller design. The sensor layers were 
given rectangular and triangular shaped electrodes, so that the charge measured is propor¬ 
tional to the tip deflection and the tip angle, respectively. The actuator layers were also 
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assumed to be covered with rectangular and triangular shaped electrodes, with the follow¬ 
ing motivation: A voltage supplied to an actuator with rectangular (or triangular) shaped 
electrodes acts in the same way on the structure as a bending moment (or force) at the tip 
of the beam. The key issue of m was to devise a stable feedback control model for that 
beam, such that it evolves asymptotically (as t —>■ oo) as a prescribed reference trajectory. 
More precisely, that controller allows to track the position and the angle of the tip mass 
at the same time. To solve the trajectory planning task, the concept of differential flatness 
(cf. [3]) was employed. Thereby, the control inputs and the beam bending deflection were 
parametrized by the flat outputs and their time derivatives. The boundary controller con¬ 
structed there has a dynamic design, thus coupling the governing PDEs of the beam with 
a system of ODEs in the feedback part. In order to render the system experimentally and 
technically realizable, it is crucial that the controller only involves boundary measurements 
up to the first spatial derivative - at the (small) price of loosing exponential stability (as we 
shall see here below). 

The goal of the present paper is first to complete the analysis of [3T], proving that this 
hybrid system is asymptotically stable but not exponentially stable. This part is an extension 
of Rao’s analysis m to dynamic controllers and inhomogeneous beams. In our second, and 
in fact main part we shall develop and analyze a dissipative finite element method (FEM) 
for the control system. 

Now we specify the problem under consideration, an inhomogeneous EBB of length L, 
clamped at the left end x = 0 , and with tip mass, moment of inertia, and boundary control 
at X = L. In the following linear system (HDl-dlSl), we actually consider the evolution of 
the trajectory error system. So, u(t, x) denotes the deviation of the actual beam deflection 
from the desired reference trajectory. Similarly, 0 i, 2 (t) denote the difference between the 
applied voltages to the electrodes of the piezoelectric layers and the ones specified by the 
feedforward controller. 


(1.1) 

fl{x)utt + {A{x)Uxx)xx 

= 0, 

0 < X < L, t > 0 

( 1 . 2 ) 

u{t, 0) 

= 0, 

t > 0, 

(1.3) 

Uxit,0) 

= 0, 

t > 0, 

(1.4) 

Juxtt{t, L) + {Auxx){t, L) 01 (t) 

= 0, 

t > 0, 

(1.5) 

Muuit, L) - {Auxx)xit, L) + 02(<) 

= 0, 

t > 0. 


Here, p, € (7^(0, L] denotes the linear mass density of the beam and A G C"^[0, L] is the 
flexural rigidity of the beam. Both functions are assumed to be strictly positive and bounded. 
M and J denote, respectively, the mass and the moment of inertia of the rigid body attached 
at X = L. Equation da states that the beam bending moment at x = L (i.e. A{L)uxx{t, L)) 
plus the bending moment of the tip body (i.e. Juxu{t, L)) is balanced by the control input 
— 01 . Similarly, describes that the total force at the free end, equal to shear force at 
the tip (i.e. —{Kuxx)x{t, L)) plus the tip mass force Mutt, cancels with the control input 
02 . 

The proposed control law has the goal to drive the error system to the zero state as 
t —>• oo. It reads: 


(Ci)t(i) = AiCi{t) +biUxt{t,L), 

(C2)t(t) = A2C2(t) + b2Ut(t, L), 

0i(t) = kiUx{t,L) + Cl ■ Ci{t) + diUxt{t,L), 

02(t) = k2u{t,L) + C 2 ■ C2{t) + d 2 Utit, L), 


with the auxiliary variables Cl, C 2 S C([0, 00 ); R") and 0 i ,02 G C'[0,oo). Moreover, A 1 , A 2 G 
l^nxn Hurwit^ matrices, fei, 62 , 01,02 G R" vectors and ki,k 2 ,di,d 2 G R. We assume 
that the coefficients ki and fc 2 are positive and that the transfer functions Gj{s) = (si — 
Aj)~^bj ■ Cj -I- dj, j = 1 , 2 satisfy 


Re{Gj{iuj)) > dj > Sj > Q Vw > 0, j = 1,2 


lA 


square matrix is called a Hurwitz matrix if all its eigenvalues have negative real parts. 
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for some constants i5i and 82 . These assumptions imply that the transfer function is strictly 


positive real, or shortly SPR (for its definition we refer to [24], [35]). Then, it follows from 
the Kalman-Yakubovic-Popov Lemma (see [53], [35]) that there exist symmetric positive 
definite matrices Pj, positive scalars e^, and vectors qj G R" such that 


^0-^3 + -^j ^3 ~ 


(1.7) 



for j = 1,2. A SPR controller is defined as a controller with SPR transfer function. One 
motivation for this controller design is the fact that, in the finite dimensional case, the 
feedback interconnection of a passive system with a SPR controller yields a stable closed- 
loop system. This principle of passivity based controller design was generalized to the 
trajectory error dynamics of the multi-layered piezoelectric cantilever in |31j . 

(fTTi) (fra constitute a coupled PDE-ODE system for the beam deflection u{x,t), the 
position of its tip u(t,L), and its slope Ux{t^L), as well as the two control variables Ci(0j 
C 2 {t)- The main mathematical difficulty of this system stems from the high order boundary 
conditions (involving both x- and t- derivatives) which makes the analytical and numerical 
treatment far from obvious. Well-posedness of this system and asymptotic stability of the 
zero state were established in [31] using semigroup theory on an equivalent first order system 
(in time), a carefully designed Lyapunov functional, and LaSalle’s invariance principle. 

In "we shall prove that this unique steady state is not exponentially stable. Let us 
compare this result to a similar system studied in m and §5.3 of |3Sj, which also consists of 
an EBB coupled to a passivity based dynamic boundary control, but without the tip mass. 
Then, that system is exponentially stable. 

As an introduction for our dissipative finite element method (FEM) in we shall now 
briefly review several numerical strategies for the EBB from the literature. In [35] the 
authors propose a conditionally stable, central difference method for both the space and 
time discretization of the EBB equation. Their system models a beam, which has a tip mass 
with moment of inertia on the free end. At the fixed end a boundary control is applied in form 
of a control torque. Due to higher order boundary conditions, fictitious nodes are needed at 
both boundaries. In [IS] the authors consider a damped, translationally cantilevered EBB, 
with one end clamped into a moving base (as a boundary control) and a tip mass with 
moment of inertia placed at the other. For their numerical treatment they considered a 
finite number of modes, thus obtaining an ODE system. In [32] the EBB with one free end 
(without tip mass, but with boundary torque control) was solved in the frequency domain: 
After Laplace transformation in time, the resulting ODEs could be solved explicitly. 

The more elaborate approaches are based on FEMs: In [5] two space-time spectral element 
methods are employed to solve a simply supported, nonlinear, modified EBB subjected 
to forced lateral vibrations but with no mass attached: There, Hermitian polynomials, 
both in space and time, lead to strict stability limitations. But a mixed discontinuous 
Galerkin formulation with Hermitian cubic polynomials in space and Lagrangian spectral 
polynomials in time yields an unconditionally stable scheme. In |13j the authors present a 
semi-discrete (using cubic splines) and fully discrete Galerkin scheme (based on the Grank- 
Nicolson method) for the strongly damped, extensible beam equation with both ends hinged. 
[3] considers a EBB with tip mass at the free end, yielding a conservative hyperbolic system. 
The authors analyze a cubic B-spline based Galerkin method (including convergence analysis 
of the spatial semi-discretization) and put special emphasis on the subsequent parameter 
identification problem. 

All these FEMs are for models without boundary control. Hence, we shall develop here a 
novel FEM for the mixed boundary control problem ([n])-([ESl). There, the damping only 
appears due to the boundary control. Hence, our main focus will be on preserving the 
correct large-time behavior (i.e. dissipativity) in the numerical scheme. Our FEM is based 
on the second order (in time) EBB equation (11.11) and special care is taken for the boundary 
coupling to the ODE. In time we shall use a Grank-Nicolson discretization, which was also the 
appropriate approach for the decay of discretized parabolic equations in [T] . We remark that 
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the modeling and discretization of boundary control systems as port-Hamiltonian systems 
also has this flavor of preserving the structure: For a general methodology on this spatial 
semi-discretization (leading to mixed finite elements) and its application to the telegrapher’s 
equations we refer to [TS] , 

The paper is organized as follows: In )J5]we first review the analytic setting from [31] for 
the EBB with boundary control. While this closed-loop system is asymptotically stable, we 
prove that it is not exponentially stable. Towards this analysis we derive the asymptotic 
behavior of the eigenvalues and eigenfunctions of the coupled system. In Sj3]we first discuss 
the weak formulation of our control system. Then we develop an unconditionally stable 
FEM (along with a Crank-Nicolson scheme in time), which dissipates an appropriate energy 
functional independently of the chosen FEM basis. We shall also derive error estimates 
(second order in space and time) of our scheme. In the numerical simulations of 2] we 
illustrate the proposed method and verify its order of convergence w.r.t. h and At. 


2. Non-exponential decay 


First we recall from m the analytical setting for (foi-dra in the framework of semi¬ 
group theory. To cope with the higher order boundary conditions (na), (na and the 
boundary terms on the r.h.s. of (na): the terms Ut{t,L), Uxtit,L) were introduced as sep¬ 
arate variables (following the spirit in earlier works [34l[^). More precisely, ip = Mv{L) 
is the vertical momentum of the tip mass and ^ = Jvx{L) its angular momentum, where 
V = Ut the velocity of the beam. Hence, we define the Hilbert space 

'H-.= {z = (u,u,Ci,C2,C,'0)"^: u e i?^(0,L),u e i^(0,L),Ci,C2 e V' e K}, 
where Hq{ 0,L) := {u S 77^(0,L)| u(0) = Ux{0) = 0}, with the inner product 

1 1 I - 1 - 

{z,z) '=2] +2 y 

+ ■^kiUx{L)ux{L) + -k2u{L)u{L) ~f PiCi -|- -(JP 2 C 2 , 
and WzWu denotes the corresponding norm. Let .4 be a linear operator with the domain 
D{A) = {z&U:u&H^{Q,L),v&Hl{Q,L)XiX2(^'^".^ = Jvx{.L),ij = Mv{L)}, 
defined by 


u 


V 

V 


77 XX 

Cl 


^iCi + j 

C2 


A2C2 + b2^ 

c 


A(Z/)'U3;3;(Z/) k\Ux(^L^ "j 

. V’ . 


(^AUxx'jxi,^') C 2 * ^2 ^‘^~M 


Now we can write our problem as a first order evolution equation: 


zt = Az 

z(0) = zo gH 


For a review of abstract boundary feedback systems in a semigroup formalism we refer 
to [5^. The following well-posedness and stability result was obtained in [31], for the 
homogeneous beam (i.e. for /r and A constant). The proof in the inhomogeneous case is 
performed analogously. Note that the contractivity of the semigroup also implies that || • ||^ 
is a Lyapunov functional for (El]). 


Theorem 1. The operator A generates a Co-semigroup of contractions on TL. For any 
zq ^ TL, has a unique mild solution z G C([0,oo);77) and z{t) Q in TL. 

But it remained an open question if this system is also exponentially stable. As a criterion 
we will use the following theorem due to Huang [53] , which was also used for controlled EBBs 
without tip mass [BES]: 
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Theorem 2. Let T(t) he a uniformly bounded Co-semigroup on a Hilbert space with infini¬ 
tesimal generator A. Then Tit) is exponentially stable if and only if 

(2.2) sup{Re(A): A € or{A)} < 0 
and 

(2.3) 

holds 

The following theorem is the main result of this section. Our proof of non-exponential 
stability of system (EH) relies on the asymptotic behavior of its eigenvalues. A related 
spectral analysis of the inhomogeneous EBB, but with a boundary control torque is given 
in m- Below we extend this study to the case when a dynamic control law is applied. 

Theorem 3. The operator A has eigenvalue pairs A„ and A„, n S N, with the following 
asymptotic behavior: 


sup ||i?(iA, A)II < oo 
AgR 


A 


n 


/(2n-l)7ry 4:hM-^pL{L)iA{L)i -1 
2h ) 2h2 


where 

(2.4) 



dw, 


and I is a real constant depending only on A, pi, and given by (12.281) . Therefore, 


sup{Re(A): A G <j{A)} = 0, 

and hence the evolution problem h2.1\) is not exponentially stable. 


Proof. We already know that the operator A has a compact resolvent (see ED). Thus, its 
spectrum a{A) consists entirely of isolated eigenvalues, at most countably many, and each 
eigenvalue has a finite algebraic multiplicity. Since A also generates an asymptotically stable 
Co-semigroup of contractions we obtain 


ReA <0, VA G a{A). 


The matrices Ai and A 2 are Hurwitz matrices and therefore only have eigenvalues with 
negative real parts. The set cr{A) fl (cr(Ai) U (j{A 2 )) C C is therefore empty or finite. Now 
we consider only such eigenvalues A of the operator A that are not eigenvalues of Ai or A 2 . 
Then z = {u, v, ((1, ( 2 , ^ D{A) is a corresponding eigenvector if and only if; 

V = Am, 

Cl = —Aux{L) {Ai — XI) ^ bi, 

C2 = -Xu{L){A2-XI)-^b2, 

and 


(2.5) {Kuxx)xx + 

(2.6) u(0) = 0, 

(2.7) u,(0) = 0, 

(2.8) A{L)uxx{T) -\- (fci — A[(Ai — XI) 61] • Cl -|- Adi + X^ J)ux{L) = 0, 

(2.9) — {Auxx)x ('^) (^2 ~ A[(A2 — XI) 62] ■ C2 -l- Xd 2 -\- X^M)u{L) = 0. 

In order to solve (I231)-(l^. we perform spatial transformations as in |21| . which convert 
(12.51) into a more convenient form. First, (12.51) is rewritten as: 

in\ I Axx , T \ 2 n 

(2.1(1) nxxxx “b ^ O.XXX “b ^ 'Oxx “b ^X U — (J. 









6 


M. MILETIC AND A. ARNOLD 


Then a space transformation is introduced, so that the coefficient appearing with u in (I2.1()|) 
becomes constant. Let u{x) = u{y), where 


( 2 . 11 ) 


y = y{x) := 


/o 


A{w) 


dw, 


with h defined as in (12.41) . Then, from (I2.6I) - (I2.10I) it follows that u satisfies: 

h'^X^u 

u(0) 


'^yyyy “1“ ^2'^yy I 


( 2 . 12 ) 


-Uy 


.(1) 


Uy(0) 

Uyy{l) + Uy{l) iPo + Kl{X)) 

■ f3lUyy{l) + /32Uy(l) + K2(A)m(1) 


with 

(2.13) 

(2.14) 


My) = 


y-xjx) 

^i{x) 


1 

2 A(a;) 


) 


- J_ /_A 

— /i2 S 16 l^A(x)J [\^A(x)J, 
/ £(£) \ 2 / 
v^y 


3 Ax(a:) 
2 A(a:) 


l^xx(x) 

A(x) 


( ifM'l ■ 

/ H(x) \ 

VA(a :)) 


( ifM'l 


and ai is a smooth function of h, and for k = 0,1,2,3. The coefficients /3o,/3i,/32 

are constants, depending on h, ^^(L), and ^^(L) for k = 0,1,2. Furthermore, we have 
introduced the following notation: 

_ 1 _ 

^i(-^) •= (x|x}) (fci - A - A/) ^ • Cl + Adi + A^ , 

(fc 2 -A((yl 2 -A/)-' 62 )-C^ + Ad^ + A^M). 

In order to solve (I2.12L we use the strategy as in Chapter 2, Section 4 of [JD]. Hence, 
to eliminate the third derivative term asUyyy, a new invertible space transformation is 
introduced: 

u{y) = 

Then (|2.12l) becomes: 

(2.15) 


(2.16) 

(2.17) 

(2.18) 

(2.19) 
where 

( 2 . 20 ) 


^VVVV ffi ^2^1/1/ + AiUy + agU + 


1 


h^X^u 

u(0) 

Uy(0) 


Uyy(l) + Uy(l) (fy + Ki(X)) + u(l) |^/34 “ -a3(l)/ti(A )J = 
-Uyyy(l) + /SgUyyil) + /dgUyil) T (/S? T K2(X)) u(l) = 


3 3 

a2(y) = a2(y) - ^asiy)^ - ^{a3)y{y), 


0 , 

0 , 

0 , 

0 , 

0 , 


and di, do are smooth functions of h, and for fc = 0,..., 4. The constant coefficients 
/33,...,/37 depend on h, ^;^(T), and ^^{L) for k = 0,...,3. Due to the invertibility of 


the above transformations, the obtained problem (I2.15I1 - (I2.19I) is equivalent to the original 
problem dLB-dlii). 

Since the eigenvalues of A come in complex conjugated pairs, and have negative real 
parts, it suffices to consider only those A in the upper-left quarter-plane, i.e. such that 
arg A G (f, 7r]. We define the unique r G C such that Re(r) > 0, and X — i'^2- It can be seen 
that argr G (0, ^]. Now, the solution to (12.151) can be approximated by the solution to the 
differential equation with the dominant terms only, i.e. Uxxxx + X^h^u = 0. More precisely, 
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we have (by adaptation of Satz 1, pp. 42 of [ 20 ]; and the last result of Lemma 12.11 is stated 
in the proof of Satz 1 ): 

Lemma 2.1. For t € (0, j], and |r| large enough, there exist linearly independent solutions 
to (12.151) . such that: 

iM = e-^-ni + / 7 (?/))- 

^ = {u:jTre‘^^^y{l + f,{y) + 0{\T\-^)), fc€{l,2,3}, 

where wi = 1 , a ;2 = i, W 3 = — 1 , W 4 = —i, and 

Jq 02 (w) dw 


fjiy) = -- 


AujjT 


+ 0{\t\ ), as |r| ^ 00 , j = 1,... ,4. 


Furthermore, the functions depend analytically on t, for j = 1,... ,4, k = 0,... ,3, 

and |r| large enough. 

Now, due to Lemma EH the solution to (I2.15|) - (I2.19D can be written as: 

w(y) = Ci'jiiy) + C'272(?/) + Csjsiy) + C'474(y), 

where the constants {Cj}'j^i are determined by the boundary conditions (12.1611 - (I2.19L 
and therefore satisfy the following linear system: 

0 = Ci7i(0)+C272(0) + C373(0) + C474(0), 

0 = Cl (7l)y(0) + 172 ( 72 )^( 0 ) + C3(73)y(0) + C'4(74)y(0), 


( 2 . 22 ) 


0 = Ei=l C'i«T-3i, 
0 = Y.i=iCim4i, 


where we define: 


1 


msi := i-f^)yyil) + i /33 + Ki(A))(7i)y(l) + i/di - -a3(l)Ki(A))7i(l), 

m^i := ~i"fi)yyyi^) + /dsi'li/yyO-) + /deili/yi^) + (Z?? + ^2(A))7i(l). 

From (I2.2ip easily follows: 

7i(0) = l + /i(0), (7^)y(0) =a;jT(l + /j(0) + O(|r|-2)), j = l,...,4, 

W 31 = e'^ + Z2 t'^)( 1 + /i(l)) + C(kP)) , 

77141 = lilsT'^ - t 3)(1 + /i(l)) + C(|r|3)) , 

77732 = {{iliT^ + Z 2 r^)(l + / 2 ( 1 )) + C>(|Tp)) , 

(2.23) 77142 = {{Ist'^ + ir^)/! + /2(1)) + Oi\T\^)) , 

mss = e"^ ii-hT^ + + /3(1)) + Cd^P)) , 

77743 = e"^ {{Ist'^ + t 3)(1 + /3(1)) + C>(|r|2))), 


77134 = e" 

77144 = e” 


( —*/lT® + l 2 T^)(l + 74(1)) + ChikP)) , 

[Ist'^ - 7r3)(l + /4(1)) + C>(|t|2)) , 


with 


_ j_ flf^Y\i^ 703(1) (lYY]\i^ := (^^Y\ 


h^A/L) \A{L) J 


hA/L) \AiL)J 


(2.24) 


= 0 . 


Ah^A/L) \AiL) J 

For u to be nontrivial, the determinant of the system (12.2211 has to vanish: 

7i(0) 72(0) 73(0) 74(0) 

(7l)y(0) (72)y(0) (73)y(0) (74)y(0) 

77731 77732 77733 77734 

77741 77742 77743 77744 

Next we shall write (12.241) in an asymptotic form when Re(r) is large: 

(2.25) i?i(7773177744 - 7774177734) + 52(7773177742 - 7774177732) + C>(|rd°) = 0, 

where 

51 := -(l + 7)[l + / 2 (l) + /3(l)] + 0 (|r|- 2 ), 

52 := (l-i)[l + /3(l) + /4(l)] + C(|r|-2). 


(2.26) 
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Noting only the terms with leading powers of r in (I2.25|l . and after division by we 

obtain 

(2.27) COST — + ^)(cosT + sin r) + 0 (|t|“^) = 0, 

4 I 3 


where 

(2.28) 


/ := 


a 2 {w) dw. 


We set fc = n — i for n S N sufficiently large and consider equation (12.271) for r in a 
neighborhood of kir. We shall apply Rouche’s Theorem (see [55], e.g.) to the equation 
(12.271) . written as 


(2.29) cosr + /(T)=0, 

where /(r) = 0(|r|“^). Consider cost on a simple closed contour K C {{n— 1)tt < Re(T) < 
mr} “around” t = kir such that |cost| > 1 on K. For n large enough, the holomorphic 
function / satisfies 1 /( 2 :) | < 1 < | cost] on K. Since t = /ctt is the only zero of cost inside 
K, Rouche’s Theorem implies that (12.2911 has also exactly one solution inside K: 


(2.30) 


T„ = fc-TT + hn 


Then, cost„ = (—I)"sin/i„. Furthermore, (12.291) implies hn = 0{n ^). To make the 
asymptotic behavior of hn more precise, we consider 


sinT„ = -(-!)”cos/i„ = -(-!)"■ +C>(n ^), 

cosT„ = (-1)"+ C>(n“^). 


Using this in (12.271) we get 

hn + + 7 ) + = 0 . 

£3 4 

Finally, this yields 

4A:7r 

and (12.301) implies 


(2.31) 



4:hPI-^fi{L)iA{L)i - I 
2/^2 


+ 0(n-i). 


Hence, condition ( 12 : 2 ]) fails and T{t) is not exponentially stable. 


□ 


In Figure [T] we show the eigenvalue pairs corresponding to the simulation example from 
im They were obtained by application of Newton’s method to the equation (12.241) . 


Remark 2.2. It can also be shown that the condition (12.311 does not hold. In particular, 
it can be shown that there is a constant C, a sequence {/in} C K diverging to + 00 , and a 
sequence {zn} C D{A) such that 


11R (f/in j A') Zn 11 T-i 

W^nWn 


C fi'fi , 


for all n large enough. 


But since the details of this calculation are rather technical we only present them in [36j . 


Remark 2.3. We shall now comment on the asymptotic behavior of the eigenfunctions of 
A. The solution to (I2.15I1 - (I2.19I1 for t = Tn has the form (see HO]): 


Un{v) 


7i(0) 72 ( 0 ) 73 ( 0 ) 74 ( 0 ) 

(7l)y(0) (72)y(0) (73)y(0) (74)y(0) 

77131 77732 7773 3 77734 

71(7/) 72(7/) 73(7/) 74(7/) 
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Figure 1. The eigenvalues A„ of the system approach the imaginary axis 
as n —>■ oo. 


up to a multiplicative constant. Using the Laplace expansion of the determinant and scaling 
the expression with has the approximate form (for n large): 


)7ry^ + + 0{n-^), 

for 0 < y < 1. Therefore, the function Un corresponding to the eigenvalue A„ has the 
following asymptotic property: 

) 7 ry^ +sin 



Unix) = e-Do”“ 3 (U<i^ -cos(^{n- 


i(y) = e 2 )’"!! _ cos ^(n - i) 7 ry^ +sin ^(n - i 


where 0 < x < L, with y = y{x) and 03 as in (12.1111 and (I2.13L 


Remark 2.4. The uncontrolled system (i.e. with Ai ^2 = 0, di ^2 = 0) is undamped and its 
operator A then has purely imaginary eigenvalues. But their asymptotic behavior is still 
like in Theorem [31 as can be verihed by the analogue of the above computation. 


3. Dissipative FEM method 

From Theorem [T] we know that the norm of the solution z(t) decreases in time. Using 
(11.71) . a straightforward calculation (for a classical solution) yields: 

^ll^llw = -^iUxt(L)^ - ^ (^(1 ■ gi + diu^t(L)J 

(3.1) — 62 UtiL)^ — - ^C2 ■ 92 + S2Ut{L)'^ 

-|c7^iCi - JCIP 2 C 2 < 0, 

where Sj = y^2(dj — 6 j), j = 1,2. Note that the r.h.s. of (13.11) only involves boundary terms 
of the beam and the control variables. Hence, ^||^||^ = 0 does not imply z = 0 (which can 
easily be verified from (ED)- 

The goal of this section is to derive a FEM for dLU-dHI) coupled to the ODE-system 
m that preserves this structural property of dissipativity. The importance of this feature 
is twofold: For long-time computations, the numerical scheme must of course be convergent 
in the classical sense (i.e. on finite time intervals) but also yield the correct large-time limit. 
Moreover, dissipativity of the scheme implies immediately unconditional stability. 
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Here we shall construct first a time-continuous and then a time-discrete FEM that both 
dissipate the norm in time. Let us briefly discuss the different options to proceed. (EH) is an 
inconvenient starting point for deriving a weak formulation due to the high boundary traces 
of u at X = L: The natural regularity of a weak solution would be u € ^([O, oo); Ho{0, L)), 
V = Ut G oo); L^(0, L)). Hence, the terms A{L)uxx{t, L), {Auxx)x{t,L) in (12.111 could 
only be incorporated by resorting to the boundary conditions El, (El- Therefore we shall 
rather start from the original second order system (fTTll-lfra. 


3.1. Weak formulation. In order to derive the weak formulation, we assume the following 
initial conditions 

(3.2a) u{0)=uoeHoiO,L), 

(3.2b) ut{0)=voeL^i0,L), 

( 3 . 2 c) Ci(o) = Ci,oeR”, 

(3.2d) C2(0) = C 2 ,o€R”. 


Moreover, let vo{L) and (vo)x(L) be given in addition to the function vq, and not as its 
trace. Multiplying El by w G ifg(0,L), integrating over [0, L], and taking into account 
the given boundary conditions we obtain; 

/'Ll pL 

/ fiUttwdx+ AuxxWxxdx + Mutt{t,L)w{L) + Juttx{t,L)wx{L) 

Jo Jo 

(3.3) +kiUx{t, L)wx{L) + k 2 u{t, L)w{L) + diUtx{t, L)wx{L) + d 2 Ut{t, L)w{L) 

-l-ci • Ci(t) Wx{L) + C 2 • C 2 (t) w{L) =0, Vic G Hq{ 0,L), t > 0. 

This identity will motivate the weak formulation. First, we define the Hilbert space 

FT :=Rx]RxL2(0,L), 


with inner product 

■■= JC(p)Ci')+MC(p)CC') + {fj,^<f,^C')L 2 , 
for ip = (^<p, '^ip, ^p), V & H. We also define the Hilbert space 

V ■={w = {wx{L),w{L),w)\ w G iFo(0, T)}, 

with the inner product 

{wi,W2)v = {A{wi)xx, {W 2 ) xx )■ 

It can be shown that V is densely embedded in H. Therefore taking H as a pivot space, we 
have the Gelfand triple 

V CH cV. 

For any fixed T > 0 we now define u = (ux{L), u{L), u) and ^i, (^2 to be the weak solution 

to (E1-(E1 and El if 

u G F^(0,r; E) niF^(0,r; H)nH^{0,T-,V'), 

Ci,C2 GiLi(0,T;R") 

and it satisfies: 


(3.4) V' < utt,w >v +a{u,w) + b{ut,w) + ei(Ci, w) -I- e 2 {C 2 ,w) = 0, 

for a.e. t G (0,T),V?h G V. The bilinear form v' < ■, ■ >y is the duality pairing between V 
and V' as a natural extension of the inner product in H. The bilinear forms a : E x E —>■ R, 
b : H X H and ei, 62 : R" x E ^ R are given by 

a{wi,W2) = {wi,W2)v+ ki{wi)x{L){w2)x{L)+ k2Wi{L)w2[L), 
b[p,v) = di{^p){^v)+d 2 {^p){^v), 

ei(Ci,w) = ci-CiWx{L), 
e2(C2,w) = C2-C,2 w{L). 
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Equation (13.411 is coupled to the ODEs 

(Cl)t(t) = + bl (^Ut(t)), 

^ ^ (C2h(t) = A2C2(t)+b2(^Ut(t)), 

with initial conditions 

(3.6a) m(0) =uo = {{uo)^{L),uoiL),uo) € V, 

(3.6b) ut{0) =vo = {{vo)x{L),vo{L),vo) € H, 

(3.6c) Ci(0)=Ci.oeR”, 

(3.6d) C2(0) =C2.o eK". 

In (|3.6all the first two components of the right hand side are the boundary traces of uq € 
Hq{0,L), but in (I3.6bp they are additionally given values. Note that in the case when 
u G H‘^{0,T; V), formulation (13.41) is equivalent to identity (13.31) . This weak formulation is 
an extension of [3] (Section 2) to the case where the beam with the tip-mass is additionally 
coupled to the first order ODE controller system. Here, we have to deal also with Ut{L) and 
Utx{L). And these additional first order boundary terms (in t), included in 5(.,.), require a 
slight generalization of the standard theory (as presented in §8 of [S^, e.g.). 

In order to give a meaning to the initial conditions (I3.6a|l , (I3.6b|l we shall use the following 
lemma (special case of Theorem 3.1 in [33] )• 

Lemma 3.1. Let X and Y be two Hilbert spaces, sueh that X is dense and eontinuously 
embedded in Y. Assume that 


u € L^{0,T-X), 
ut e L^{Q,T-Y). 


Then 

ueC{[0,T]-,[X,Y].]), 

after, possibly, a modification on a set of measure zero. Here, the definition of intermediate 
spaces as given in Ea, § 2 .1, was assumed. 

Theorem 4. (a) The weak formulation CT - (ES) has a unique solution {it, (^ 1 ,( 2 )■ 

(b) The weak solution has the additional regularity 


(3.7a) 

(3.7b) 

(3.7c) 

(3.7d) 


ueL^{ 0 ,T-,V), Ut&L^{ 0 ,T-,H), 

Ci,C 2 €C([ 0 ,T];M"), 

ueC{[ 0 ,T];[V,H].), 

Ut€C{[ 0 ,T]-,[V,H]\). 


Furthermore, even stronger continuity for the weak solution can be shown: 


Theorem 5. After, possibly, a modification on a set of measure zero, the weak solution u 
of lETP-lElP satisfies 

u G C([ 0 ,T];E), 

Ut G C{[0,T]-H). 

The proofs of Theorem |3] and [5] are given in Appendix A. 


3 . 2 . Semi-discrete scheme: space discretization. Now let Wh C Hq{0,L) be a finite 
dimensional space. Its elements are globally C^[0,L], due to a Sobolev embedding. For 
some fixed basis Wj,j = 1,...,N the Galerkin approximation of ( 13 . 41 ) reads: Find uu G 
C'^{[0,oo),Wh), i.e. Uh = {{uh)x{L),Uh{L),Uh) G (^^([O, 00), E), and C1.2 G (^^([O, 00), R”) 
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with 

Jo dx “t" Jq XX dx 

+M{uh)tt(t, L)wj{L) + J{uh)xtt {t,L){wj)x{L) 

(3.8) +ki{uh)x{t,L){wj)x{L) + k2Uh{t,L)'Wj{L) 

+di{uh)xt{t, L){wj)x{L) + d 2 {uh)t{t:L){wj){L) 

+ Cl ■ Ql{t) {Wj)x{L) + C2 • C2(<) Wji-L) = O’ j = 1, • ■ ■ i > 0, 
coupled to the analogue of (EH): 

/oq'i (Ci)t(O = AiCi(t)+ 6 i(u,,)^t(t,i), 

^ (C2)t(<) = A2C2(i)+?>2(zi.)*(t,2^), 

and the initial conditions 


M/i(0, .) 
(u/t)i(0,.) 

Ci(0) 

0(0) 


M?i,0 s 

S W/t, 

Ci.o e K", 
C2,0 6R". 


(13.81) is a second order ODE-system in time. Expanding its solution in the chosen basis, i.e. 

N 

Uh[t, x) = ^ Ui{t)wi{x), 

i=l 


and denoting its coefficients by the vector 

U = [ C/i U 2 ... Un]'^ 

yields the equivalent vector equation; 


(3.10) 


AUtt + lUt + KU + C{t) = 0. 


Its coefficient matrices are defined as 

Aij := / ^WiWj dx + Mwi{L)wj{L) + J{wi)x{L){wj)x{L), 

Jo 

Bjj := dl{w^)xiL){wj)xiL) + d2W^{L)wj{L), 

A{wi)xx{wj)xx dx + ki[wi)x{L){wj)x{L) + k 2 W^{L)wj{L), 

ij = 



and the vector C has the entries 


Cj{t) = Cl • Ci(i) (wj)x{L) + C2 • C, 2 {t) Wj{L), j = 1,... ,iV. 

The matrix IK is symmetric positive definite, since we assumed fci _2 > 0. Since also A is 
symmetric positive definite, one sees very easily that the IVP corresponding to the coupled 
problem ()3.10ll . (13.91) is uniquely solvable. 

For a final specihcation of the FEM we need to choose an appropriate discrete space. 
Only for notational simplicity, we shall assume a uniform distribution of nodes on [0, L]: 

Xm = rnh, m € {0,1,..., P}, 

where h = p. A standard choice for the discrete space Wh is a space of piecewise cubic 
polynomials with both displacement and slope continuity across element boundaries, also 
called Hermitian cubic polynomials (see [23], [B], e.g.). They have been employed not only 
for the Euler-Bernoulli beam, but also Timoshenko beams (cf. [Hj). To define a basis for 
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Wh (Hermite cubic basis, see e.g. [IS]), we associate two piecewise cubic functions with each 
node Xrm m>l satisfying: 

W2m-i(a;fc) = I q’ ^m-i(a;fe) = 0, 

w' 2 m{xk) = ^^Q W 2 r„(a;fc)= 0 , 

for all k = 0,..., P. Hence, the nodal values of a function and of its derivative are the 
associated degrees of freedom. Due to the boundary conditions at x = 0 in Wh C Hq, the 
basis set does not include the functions W-i and wq associated to the node Xq = 0. Thus, 
N = 2P. For the coupling to the control variables we shall need the boundary values of 
Uh- The above basis yields the simple relations Uh{t,L) = 17jv-i(i), {uh)xit,L) = UN{t). 
Compact support of the basis functions leads to a sparse structure of the matrices 

A, B, and K: A and K are tridiagonal, B is diagonal with only two non-zero elements 
BAr_i,Ar_i = c? 2 , BAr,Ar = di. And the vector C has all zero entries except for Civ_i = C 2 ■ C 2 , 
Cat = Cl • Cl- 

Next, we shall show that the semi-discrete solution Uh{t) decreases in time. As an analogue 
of the norm || 2 ;(t)||'H from [21 we first define the following time dependent functional for a 
trajectory u G C'^{[0 ,oo)]Hq{0 ,L)) and Ci ,2 G 00 ); K"): 

1 M J 

A(t;u, ( 1 ,( 2 ) '=2 7 {^Ua:x{t,xf'+ ^iutit.xf') dx++—Uxtit,LY 

(3.11) +-^Ux{t,LY + -^u{t,L)‘^ + -Cj (t)PiCi(t) + {'t)P2C2{t). 

For a classical solution of (EH) in D{A) we have P(t; u, Ci, C 2 ) = lk(0llw- 

Theorem 6. Let Uh G cx)); Pq (0, P)) and Ci ,2 G cx));M") solve (13.81) . (13.91) . 

Then it holds for t > 0: 

^E{t]UhXiX2) = -^CiPiCi-^{Ci-qi + Si{uh)xtiL)^ -Si{uh)xt{Lf 

~~^^2 P 2 C 2 - 2 (c 2 ■ <?2 + hiuh)t{L)"j - d 2 {uh)t{LY < 0. 
Proof. In the following computation we use (13. 8|) with the test function Wh = {uh)t- 
d f ^ 

j Cl; C 2 ) — / ^{.'^h^xx{ph)xxt dx -\- I dx 

+M{uh)t{L){uh)tt{L) + J{uh)tx{L){uh)ttx{L) 

+ki{uh)x{L){uh)xt{L) + k2iuh){L){uh)t{L) 

+ c7 -Pl(Cl)t + <^ 2 '^ ^2(C2)t 

= -di{uh)xt{LY - d2{uh)t{Lf 

— Cl • fl{Uh)xt{L) — C2 ■ f2{uh)t{L) + c7-Pl(Cl)t + C2^ P2iC2)t, 
and the result follows with (13.91) and (11.71) . □ 

In the undamped case (i.e. Aj = 0,dj = 0) the energy E is clearly preserved in the 
semi-discrete system. Furthermore, it has been shown in the proof of Theorem [5] that the 
energy functional for the weak solution u, Ci i C 2 of E3) - EH) has an analogous dissipative 
property, cf. (15.131) . 


3.3. Error estimates: semi-discrete scheme. Since using cubic polynomials for the 
space approximation, we shall obtain accuracy of order two in space (in P^(0, L)). Thereby, 
the common method for obtaining error estimates (cf. [T3]) will be adjusted to the problem 
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at hand. With u we denote the nodal projection of the weak solution u to Wh, defined in 
terms of Hermite polynomials: 

p p 

U{t,x) = +'^Ua;{t,Xm)w 2 m{x). 

m—1 m—1 

Assuming that 

ueC([0,T];ff4(0,A)), 

(3.12) utGL^iO,T;HliO,L)), 

Utt€L^i0,T-,H^{0,L)), 
it can be seen (e.g. in [7], [13]) that a.e. in t\ 

\\u — u\\h^(0,L) < Ch'^\\u\\H‘>-{0,L), 

(3.13) \\ut — Ut\\H^(o,L) < 

\\utt — Utt\\L^{0,L) < Ch'^\\utt\\H^{0,L)- 

We define the error of the semi-discrete solution {uh, ( 1 ,( 2 ) as et '■= Uh — u £ Wh and 
Cf := Ci ~ Cij * = li 2. Then using (I3.8I1 - (I3.9I1 we obtain 


Jq fj-ieh)ttw dx + Jq A{eh)xxWa:xdx 
+M{eh)tt{t,L)w{L) + J{eh)xtt{t,L)wx{L) 

+ki{eh)x{t,L)wxiL) + k 2 eh(t, L)w{L) 

+di{eh)xt{t, L)wx{L) + d2{eh)t{t,L)w{L) 

-fci • Cr(i) Wx{L) + C 2 ■ C|(i) w{L) 

= /(, n{utt - Utt)w dx + /g A{uxx - Uxx)wxx dx, Vic e Wt, t > 0, 

coupled to: 


(Ci)t(i) — AiQ{t) + bi{eh)xt{t, L), 

(Cl)t(i) = A2Ci{t) + b2{eh)t{t,L). 

Using w = {eh)t and proceeding as in the proof of Theorem[6]we obtain 

< Jq - utt)ieh)tdx + A{uxx - Uxx)ieh)txxdx, 

for a.e. t £ [0,r]. Integrating (13.1411 in time, and performing partial integration, we get 

E{t;ehX!XI) < E{0-,eh{0)XmXm) 

+ 2/g/g ix{uu{s,x)-utt{s,x)){eh)t{s,x)dxds 

(3.15) 2 Jg A{v,xx{d^x) Uxxid ^ xYjich^xxidi dx 

“f 2 Jq A(zi3,a;(0, x) Uxx{^ix')'){c}i)xx{^iX^ dx 

^ So So x) Utxxi^j x')')(^ 6 fi^xxi^: X^ dx ds. 

Applying Chauchy-Schwarz to (13.151) yields: 

(3.16) 

Eit-eu.CiXl) < i?(0;e,,(0),Cf(0),CI(0)) 


Eflmax ^tt|lL2(0,T;L2(0,L)) /o 

-fA 

77 iax ( 8 || Uxx(t, ■) '^xx{t, .)\\L^(o,L) + 8 ||(£?i)xx(^, OIIls^o,^) 

+ 8||Mxx(0, .) — Uxx{Q, ■)\\l^(0^L) 8 II ■) IIl^(0,L) 

+ \\ut — '“t|lL2(0,T;H2(o,L)) + So II (£?i)2:a:(s, •)IIl2(0,L) J 

where fJ,max = Tx>-ax.x^[o,L] A^ax = xa.a,yix^[o,L] A{x). Next, we use (13.131) to obtain: 

|U(t; ehXl, Cl) < e/.(0), Cf (0), Cl(0)) + 2 f* E(s; Cl, Cl) ds 

+ Ch‘^ (l|u||p([g 7.].^4(o^i)) + lki|lL2(o.r;iL4(o,L)) + ll^«lll,2(o_T;iL2(o,L))) ' 


(3.17) 
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Gronwall inequality applied to (13.1711 gives: 


(3.18) 


E{t- Ch, Cf, Cl) < C (£;(0; 6^,(0), Cle(O), C2e(0)) 

+ ( ll“llc'([0,T];ff<‘(0,L)) + ll“‘lli2(0,T;ff*(0,L)) 


WuttW 


2 

L2(o,T;ff2(0,L)) 



Finally, we have: 

Theorem 7. Assuming h3.12\] . the following error estimate of the semidiscrete solution 
holds: 

E{t- Uh - u, Cl - Cl, C2 - C2)^ < C (e( 0; eh(0), Cle(O), C2e(0))5 

(3.19) + (lluttllL^(0.T;H^(0,L)) + \\ut\\L^{0,T-H^{0,L)) + ||'u||c([0,T];iI4(o,L)))^ , 

0<t<T. 


Proof. The result follows from (13.131) . (I3.18p . and the triangle inequality. 


□ 


3.4. Fully discrete scheme: time discretization. For the numerical solution to the 
ODE (13.101) we first write it as a first order system and then use the Crank-Nicolson scheme, 
which is crucial for the dissipativity of the scheme. To this end we introduce Vh '■= {uh)t, 
and V := Ut = [ Fl V2 ... Vn is its representation in the basis {wj}. The solution of the 
system (13.811 . (13.911 is then the vector Zh = [uh Vh Ci C 2 ]^- In contrast to (JU here we do not 
have to include the boundary traces Vh{L), (vh)x[L)\ In the finite dimensional case both Uh 
and Vh are in Hq{0, L). In analogy to [J2l the natural norm of Zh = Zh{f) is defined as 

(3.20) ||z,f := -J^A{uh)lxdx + -j^^,vldx+-vl{L) + -{vh)l{L) 

+ Y(n/i)x(-I') + + ^Ci^-^iCi + 2CJ-P2C2- 

Let At denote the time step and 


tn = nAtjVn G {0,1,..., S}, 

is the discretization of the time interval [0,T], T = S'At. For the solution of the fully 
discrete scheme at t = tn, we shall use the notation z" = [u" u" C" C 2 ]^- ^nd U", V" are the 
basis representations (in {wj}A_^) of u" and v”, respectively. Furthermore, let the vector 
C" be defined by: 


(C"), := Cl • Cr {u>ML) + C 2 ■ C2” w,{L), j = 1,..., A. 
The Crank-Nicolson scheme for (13.1011 . (13.911 then reads: 


(3.21) 


(3.22) 

(3.23) 

(3.24) 


U"+i _ u" 


At 


AV”+i - AV” 

At 


^n+1 

Si 

cr 

At 


^n+1 

S2 

Cl 


At 


i(V"+i+V”), 

-i(KU”+^ -f KU”) - i(]BV"+^ -f BV”) 

-i(C”+^ -f C”), 

<+^(L)+<(L) 

^ 2 ^ 2 


In the chosen basis {rcj}, the last term of (13.231) . (13.241) reads + Vf^') /2 and (V^^i + 

respectively. Next, we show that this scheme dissipates the norm. The somewhat lengthy 
proof is deferred to the Appendix B. 
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Theorem 8 . For n S Nq it holds for the norm from (13.201) ; 




C-+1 + cr 


2 191- 


<5i 


At 


<+i(L)-<(L) 


2 ' At 


At 


+ 2U2- 


+ S2 


u^+^{L) -u^{L) 
At 


//*"+! I /-n\ I ^n+l I An , _l_ /‘"I I /‘W+l i An 

, ei (Cl +Clj n Cl +Cl , £2 (C2 +C2j n C2 + C2 

+ -1--r t;--^2-- 


This decay of the norm is consistent (as At 0) with the decay (13.11) for the continuous 
case, and with Theorem [51 For the uncontrolled beam (i.e. 0i = 02 = 0), Theorem [5] 
shows that ||z"|| is constant in n. This motivates our choice of the Crank-Nicolson time 
discretization. 


Remark 3.2. Note that the scheme (I3.21|) - (|3.24ll and the norm dissipation property from 
Theorem IS] were written independently of the basis {wj}. Hence, this decay property applies 
to any choice of the subspace Wh C Hq{0,L). And the same remark applies to Theorem [6l 

3.5. Error estimates: Fully discrete scheme. In this subsection we shall need to assume 
additional regularity of the weak solutions u, Ci and C 2 , in order to estimate the error of the 
fully discrete case: Suppose that u € Hq{0, L)) and Ci:C 2 G iJ^(0, T; R"). Let us 

define u € Wh to be the projection of the weak solution u, such that 

a{u{t),Wh) = a{u{t),Wh), 

'it € [0,T]. One easily verifies that it holds: u G H'^{0,T; Hq{0, L)), since the projection 
u I— > it is bounded in F[q(0,L). Furthermore, let u® := u — u denote the error of the 
projection. Assuming u G H^{0,T; IIq{0, L)), we obtain the error estimates for ii (cf. [IS]): 

ll'a'^l|ff2(o,L) < C'ii^||u||//4(o,L), 

(3-25) ilwi j|_f/2(o_L) < C'ii^||ut||//4(o^i), 

l|ntt||ff2(o.L) < Chf\\utt\\H^[oA)- 

Let z{tn) = [u{tn) v{tn) Ci(in) C 2 (ira)]^ and z" = [u" u" Cff QY denote the solution of the 
system and the solution of the fully discrete scheme at time i = i„, respectively. Then we 
define the error by 

e" := u” - u(tn), 

:= v"‘-ut{tn), 

■■= C-C^{tn), *= 1 , 2 , 

and := [e’" Qi (” 2 ]^. for every n G 0,1,..., S'. 

We now give the second order error estimate (both in space and time) of the fully discrete 
scheme. The proof is deferred to Appendix B. 

Theorem 9. Assuming u G H'^{0,T; Hq{ 0, L)) Ci F['^{0,T-, Hq{ 0, L)) and 
Ci,C 2 € i?^(0, T;M"), the following estimate holds: 

Ik” “ ^(fra)ll ^ C” [||Zg|| + il^||'u||_fL2(0_T;ff*(0.L)) + (^0^ ( 11**** || 'i (O.L)) 

+ ll***t||ff2(o,T;ff2(o__L)) + ||(Cl)*t||L/i(0.T;B") + || (C2)t* ||ff * (O.T;R'*)) ] ■ 


4. Numerical Simulation 


In this chapter we verify the dissipativity of our numerical scheme for an example with the 
following coefficients: g, = A = L = l, M = J = 0.1, ki = k2 = 0.01, and di = d2 = 0.02. 
We take n = 10 as the dimension of controller variables. Thereby, Ai = A 2 = —I G 
where I is the identity matrix, and 61 = 1)2 = ci = C 2 = [1 1 ... 1]^ G We take the 
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Figure 2. Damped vibration of the beam: deflection u{t,x) 


time step At = 0.01 and the spatial discretization step h = 0.01. Figure [2] shows the damped 
oscillations of the beam and its convergence to the steady state m = 0 on the time interval 
[0,50]. Figure [3] illustrates the (slower then exponential) energy dissipation of the coupled 
control system. Finally, we perform simulations for different time and space discretization 



Figure 3. Dissipativity of the norm (or “energy”): ||z(<)||-h 


steps to verify the order of convergence (o.o.c.) proved in UH In Tabled] we list the Z^-error 
norms of Zg- In the left table we see the o.o.c. results for fixed At = 0.01 and varying 


Table 1. Experimental convergence rates 


At 

h 

ll-^e 11/2 

O.O.C. 

10-2 

1 

f 

! 

1.75*10-2 

— 

10-2 

5.5*10-3 

1.67 

10-2 

7.92 * 10-"^ 

2.80 

1 

o 

1 — 1 

1 

¥ 

64 

1.39* 10-4 

2.51 

10-2 

3.38* 10-3 

2.04 

10-2 

1 

128 

8.24* 10-3 

2.04 


At 

h 


O.O.C. 

6.4* 10-3 

1 

50 

2.58* 10-3 

— 

3.2*10-3 

1 

50 

6.87* lO-"^ 

1.91 

1.6* 10-3 

1 

50 

1.73* lO-"^ 

1.99 

8* 10-^ 

1 

50 

4.27* 10-3 

2.02 

4* 10-^ 

1 

50 

1.02* 10-3 

2.07 

2* 10-^ 

1 

50 

2.03* 10-® 

2.32 
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space discretization step h on the time interval [0,1]. In the right table the o.o.c. results for 
different At but h = 1/50 fixed, on the time interval [0,0.00041], are presented. 

5. Appendix A 

The following proof is an adaption of the proof of Theorem 8.1 in m, for the system 
studied here. It is included for the sake of completeness. 

Proof of Theorem^ (a)-existence: Let {'Wk}^i be a sequence of functions that is an or¬ 
thonormal basis for H, and an orthogonal basis for V. We introduce Wm '■= span{wi ,..., thm}, Vm G 
N. Furthermore, let sequences UmO,VmO G Wm be given so that 

un in V, 

(5.1) 


UmO 

VmO Vo in H. 

For a fixed m G N we consider the Galerkin approximation 


Vm(t) — j Umi^L^ j Um^ — ^ ^ ; 


fc=l 

with d^{t) G K, which solves the formulation (13.81) for all w G Wm'- 

(5.2) {{Um)tt,w)H + a{Um,w) -f b{{Um)t,w) -f ei(Cl.m,Mi) + e2((2,m,w) = 0, 

and Ci,m,C 2 ,m solve the ODE system 

/c o'! (Cl,m)t( 

^ (C2.,n)t( 

with the initial conditions 


— ^lCl,'m 

it) 

+ bi^ 


= ^2C2,m 

it) 

+ ^ 


Vm ( 0 ) 

= 

^mO 

•) 

(Um)t(O) 

= 



Cl,m( 0 ) 

= 

Co,i, 


C 2 .m( 0 ) 

= 

Co,2- 



This problem is a linear system of second order differential equations, with a unique 
solution satisfying Um G C^([0,T];E) and Ci,m,C 2 ,m G (^^([O,Tj;R”). Next, we define an 
energy functional, analogous to (I3.11|l . for the trajectory {uXiX^)' 


Eit]uXlX2) ■■= ^l|M(i)lly + Y(^M(t))2 -b y(2M(t))2 


^2 ,2 ■ 


\\Mt)\\ 


+ 2^1^ (^)^lCl(^) + 2 ^ 2 '' {t)P2C2{t) 


= II (m, Ut, Cl, C2, Jutx{J), Mut{L))\\-H. 

Taking w = (um)t in (15.21) and using the smoothness of Um, Ci.m, C2,m, a straightforward 
calculation yields 

— 5i{}{Um)t)^ ~ 2 (Cl,m ■ 9l + ^l(^(Mm)t)^ 
-52{^{Um)tf- 


dt 


A(t, 1/77.1, Cl,m, C 2 ,m) — 


2 (C2,m ■ 92 + d2C{Vm)t 
e 

,m 

F(t^ Urm Cl,in; C 2 ,m) ^ 0 


--^(Cl,m)^^’lCl,m- ^iC 2 ,m)^ P2C2, 


(5.4) 

which is analogous to (13.11) for the continuous solution. Hence 

.^(t, l/m, Cl,?^^, C 2 ,m) A A( 0 , 1/771, Co, 1 , Co, 2 ), t ^ 0 , 


which implies 
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(5.5) 


{Cs.mlmGN 


is bounded in C([0,T];F), 
is bounded in C{[0,T\,H), 
are bounded in (^([O, T]; R"). 


Due to these boundedness results, it holds Vu) € V: 

la(u^(t),w) + b((Um)t(i),w) + ei(Ci^rn(t), w) + e 2 (( 2 ,ni(t),w)l < Di\\w\\v, 

a.e. on (0, T), with some constant Di > 0 which does not depend on m. Now, let m G N be 
fixed. Furthermore, let w gV, and w = wi + W 2 , such that wi € Wm and W 2 orthogonal to 
Wm in H. Then we obtain from (15.21) : 

= -a{Um,Wi) - b{{Um)t,Wi) - ei(Cl,m,W>l) - e2(C2,m, Wl) 

< Di\\wi\\v < Di\\w\\v. 

This implies that also {um)tt is bounded in L^(0,T;D'). Furthermore, from (15.31) it trivially 
follows that {(Ci,m)t}mGN and ({C 2 ,m)t}mGN are also bounded in L^(0,T;M"). 

According to the Eberlein-Smuljan Theorem, there exist subsequences {um, }zgn, {Ci.mi }zeNi 
{C 2 ,m,}zGN, and u G L'^{0,T]V), with ut G L'^{0,T\H), uu G L'^{0,T]V'), and Ci,C 2 S 
(0, T; R”) such that 


(5.6) 


{Umi } 
{{Umi)t} 
{{Umi)tt} 

{Cl,mi} “ 

ICz.mi } “ 
{(Cl,mi)t} “ 

{(C2,mi)t} 


‘ u in L2(0,T;D), 

-Ut in L^iO,T-,H), 
-utt in T2(0,r; V), 
Cl inL2(0,T;R"), 

C 2 in L2(0,r;R"), 
(Ci)t in L2 (o,T;R"), 
■(C 2 )t in L2(o,r;R"). 


Therefore, passing to the limit in (15.21) and (15.31) . we see that u and C 11 C 2 solve (13.41) and 
(1^ . 

(^6}“additional regularity: From C 11 C 2 S Ff^(0,r;R"') follows the continuity of the con¬ 
troller functions, i.e. (I3.7bl) . It is easily seen from the construction of the weak solution and 
(15.51) that u satisfies (|3.7ap . (I3.7cl) follows immediately due to Lemma ITTI after, possibly, 
a modification on a set of measure zero. (I3.7dl) follows from Lemma 13.11 and the ’Duality 
Theorem’ (see [33|, Chapter 6.2, pp. 29) which states: for all 6 G (0,1), it holds 

[x,r]'=[r',x']i_,. 

(a)-initial conditions, uniqueness: It remains to show that u, Cij nnd C 2 satisfy the initial 
conditions. For this purpose, we integrate by parts (in time) in (13.41) . with w G C^OO, T]; V) 
such that u)(T) = 0 and Wt{T) = 0: 


(5.7) 


/ [{u,Wtt)H + a{u,w) + b{ut,w) -I- ei(Ci, w) + e 2 (C 2 ,li)] dr = 

Jo 

-{u{0),wt{0))H + V' < Ut{0),w{0) >v- 
Similarly, for a fixed m it follows from (15.21) : 

/ [{Um, Wtt)H + a{um, w) + b{{um)t, w) + ei(Cim, w) -|- e2(C2m, w)] dr = 

Jo 


(5.8) 


-{UmO,Wt{0))H + {VmO,w{0))H- 
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Due to (EH) and (EH). passing to the limit in (15.811 along the convergent subsequence {umi} 
gives 

/ [{u,Wtt)H + a{u,w) + b(ut,w) + ei(Ci,w) + e 2 {C, 2 ,w)] dr = 


(5.9) 


-{uo,WtiO))H + {vo,w{0))h- 


Comparing EH with (15.91) . implies {t(0) = uq and 'Ut(O) = vq. Analogously we obtain 
Ci(0) = Co,i and C 2 ( 0 ) = Co, 2 - 

In order to show uniqueness, let {u, Ci, C 2 ) be a solution to (15.4|) and (13.51) with zero initial 
conditions. Let s G (0, T) be fixed, and set 


and 


^i(i) := f (iir)dT, 

Jo 


for * = 1,2. Integrating (13.511 over (0,t) yields with (11.711 


(5.10) 


+ {di - Si){’'u{t))'^ + Zi{t) ■ Ci(*M(t)), 


for 0 < t < T, i = 1, 2. Integrating (13.41) with w = U over [0,T], and performing partial 
integration in time, yields 

[ {ut{T),u{T))H - a{Ut{T),U(T)) + b(u{T),u{T)) dr 


2 

(5.11) / Zi{T) ■ Ci{''u{T)) dr = 0. 


From (15.101) and (I5.11|l follows 


Jo -^a{UiT),U{T))+ ^Y^^ZJ{t)P,Z,{t)^ dr 

2 


2=1 ' 


(5*(*n(r))2 + {t)P,Z,{t) + ^{qi ■ Z,{t) + S^{^u){t))'^ ) dr. 


Therefore, 


^ll«(s)||?f + ia(y(0),{7(0))+^iz7(s)P,Z.(s) < 0. 

i=l ^ 

The matrices Pj,j = 1, 2 are positive definite, and the bilinear form a(.,.) is coercive. Hence 
u{s) = 0, 17(0) = 0, and Zi{s) = 0. Since s G (0,T) was arbitrary, m = 0, Ci = 0, z = 1,2 


follows. 


□ 


Before the proof of the continuity in time of the weak solution, a definition and a lemma 
will be stated. 

Definition 5.1. Let F be a Banach space. Then 

C„([0,T];y) := {zn G L“(0, T; F): V/G F' 

t i-A {f,w{t)) is continuous on [0,T]}. 

denotes the space of weakly continuous functions with values in F. 

The following Lemma was stated and proven in [33] (Chapter 8.4, pp. 275). 
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Lemma 5.2. Let X, Y be Banach spaces, X G Y with eontinuous injection, X reflexive. 
Then 

L^fl),T-X)GC^fl),T- Y) = C^{{),T-X). 

Proof of Theorem [5l This proof is an adaption of standard strategies to the situation 
at hand (cf. §8.4 in [53] and §2.4 in jUj). Using Lemma [5.21 with X = V, Y = H, we 
conclude from (I8.7all . (I8.7cl) that u S Cu,([0, T]; U). Similarly, (|3.7all and (I3.7d|) imply 
UteC^i[0,T]-,H). 

Next, we take the scalar cut-off function Oj € such that it equals one on some 

interval I CC [0, T], and zero on IR\ [0, T], Then the functions Oju : K U and O/Ci, O/C 2 : 
R —>■ R” are compactly supported. Let 77 *^ : R —?► R be a standard mollifier in time. Then we 
define 


r := r]^ *Oiu€C^{R,V), 

Cl ■■= 77^*0/Ci ecr(K,K"), 

C 2 ^ := 77^*07(2 €C'“(R,R"). 

Now Cl and Cl converge uniformly on I to Ci and C 2 , respectively. Moreover, u'^ converges 
to u in V, and ul to Ut in H a.e. on I. Then, .E(t; m'^, Ci, Cl) converges to E{t-,u, ( 1 X 2 ) a.e. 
on I as well. Since u'^jCijCI are smooth, a straightforward calculation on I yields 

(5.12) ^U(t;7i%C|,CI) = F{t-,u^,C,C), 

with F defined in (15.41) . Passing to the limit in (|5.12l) as e ^ 0 

(5.13) ^E(t-,u, ( 1 X 2 ) = F{t-,uXiX2) 

at 

holds in the sense of distributions on I. Since I was arbitrary, (15.1311 holds on all compact 
subintervals of (0, T). Now 1 1 —>■ E{t; u, Ci, C 2 ) is an integral of an L^-function (note that the 
input functions of F satisfy: lit G L^(0,T)), so it is absolutely continuous. 

For a fixed t, let lim„_>+oo X = t and let the sequence Xn be defined by 

Xn ■■= - Htn)\\v + ^WMt) - MXWh 

+ y(^w(0 - + y (^u(i) - ^u(tn)f 

+ lXlit) - Citn)V PlXlit) - Cl (in)) 

+ 2 (‘^ 2 (t) - C2(tn)VP2X2it) - Xitn))- 

Then 

Xn = F(t; u, ( 1 X 2 ) + F(X;u, ( 1 X 2 ) - (u(t),u(tn))v - (ut(t),Ut(tn))H 

-ki ^u{tCu{tn) - k2‘^u{t)‘^u{tn) - CXV PlC{X) “ C 2 (i)^f2C2 (in) • 

Due to the t-continuity of the energy function, weak continuity of u,ut, and continuity of 
Ci,C 2 , it follows 

lim Xn = 0. 

n— )-+oo 

Finally, it follows that 

lim ||7it(t) - 7it(t„)||^ = 0, 

n—^oo 

lim \\u{t) - uiflXWv = 0) 


which proves the theorem. 


□ 
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6. Appendix B 

Proof of Theorem O First we obtain from (13.2111 and (13.2211 (written in the style of (13.311 i: 


( 6 . 1 ) 


At 


( 6 . 2 ) 


rL ^n+l _ 


At 


■Wh dx - 


rL n+1 I n 
^ ^XX ' ^xx 


it^h^xx dx 


, w"+i(A)-i;"(L) v^+\L) -v^{L) ^ ^ 

+ M -- Wh{L) + J - —- (wh)x{.L) 


+fcl 


At 

<+i(L)+<(L) 


At 

r ^ rr, , u'^+\L) + U^{L) . , 

{wh)x{L) + k 2 -r- -Wh{L) 


. <+i(L) + <(L), , , v^+\L) + v^{L) , . 

-\-di -^- (wh)x{L) + d2 - - - Wh{L) 

+Ci • — — —^{wh)x{L) + C2 • —— - -Wh{L) = 0, ywh G Wh- 


Next we multiply (EUl by — u”), and integrate over [0,i] to obtain 


i dx= r\!^—^iv-+^-vndx, 


and Wh = in (16.211 : 


Huxt^fdx = -- / fi 


L ^n +1 _ 




At 


dx 


, u"+i(A)-z;"(i) 
-M --AAu 


"+i(L)-J 




At ^ ' At 

n+^( T\ ,,n( T\ .,,"+11 


ity' *(L) 

-t. <« (L) - (L) 

-lii- 




-Cl 


>n+i j_ /-n+l . 

”*”^1 n+l/7-\ Ss2 ^2 n+l/rx 

-^ [L) - C 2 - - - u ^ [L). 


We next set Wh = u” in (16.211 : 


HKx)'^ dx =-^ H 


L ^n +1 _ 


At 


' dx 


-M- - ^ ^ ^ u"(L) - J " ^ ^ <(L) 


-fci 


At 

<+>(L) + <(L) 


At 




-Cl • 


(-+1 + 


/•^+1 _[_ /-n 

S2 S2 


U-{L)-C2- ^^ 7"^ u"(A). 


2 


2 
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This yields for the norm of the time-discrete solution, as defined in (13.201) : 


,ri+l||2 


II2 




v^+^(L)^ - v^{Ly 
2 

- v:{Ly‘ 


At 


+ 


2 

k 2 


(- «+i(L) + <(L)) «+i(L) - <(L)) + u:+\Lf - <(L)2) 
(- {u^+\L) + u^{L)) (u'-^+^L) - u'^iL)) + u^+\Lf - u'^ihf) 


^«+i(L) + u:(L))«+i(L) - <(L)) 

^(u"+i(L) + u"(L))K+i(L)-u'‘(L)) 

ici • (cr+'+cr)«+'(i) - <(i)) + ^(cr')^Picr' - ^(cD^^^icr 
ic2 • (Cr^ + C2")(^”+'(i^) - + \{C2^^VP2Q^^ - \{QVP2Q. 


For the first six lines we use (I3.2ip . and for the rest Cj = Pjbj -I- qjSj (cf. (11.71) ') to obtain: 


n+l||2 _ 


|2"lP- 


««(l) - <(i))“ - ||(«”'(r) - “"(i))" 


(cr+^ + cr) 

2 

(C2+ C2") 


■(Fi&i+gi<5i)«+'(L)-<(L)) 


■(P2&2+g2<52)(u"+nT)-u”(T)) 


+ licrYPicr^ - ^(cD'^Picr+^(C2+')^^2C2”+' - lic^VP2c^. 


(6.3) 


For the second and the third line of (16.31) we now use (13.211) . (13.231) . and (13.241) from the 
Crank-Nicholson scheme: 


^n+l||2 


^ {<^\L) - <{L)f - - u-{L)f 

(cr+' + cyf 


(cr+'+ cr) 


cr + cr+' 


Pi cr' - cr - At A: 


q^ 6 ,iur\L)-u:iL)) 


(cr+^ + cr) ^^ cr+^ + cr 

(cr^^ + cr) 


q252ir+\L)-riL)) 

+ i(cr+')^^icr+' - i(cr)^Picr + ^(C2”+')^^2cr+' - l{cr^P2C2- 
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Since Pj , j = 1, 2 are symmetric matrices, this yields 


di 




(cr+' + cr) 






+ C2) 


■q2S2iu^+\L)-u^{L)), 


which is the claimed result (by using (ll.7|) V 


□ 


Proof of Theorem Let k € {0,1,...,5'} be arbitrary. Taylor’s Theorem yields Vx G 

[0,L]: 


(6.4) 

where 


u(tk+i,x) - u{tk,x) _ Ut{tk+l,x) + Ut{tk,x) 
At ~ 2 


+ Atr'=(x), 


T,^ix) = 


^ (4+1 -dt + / (4 - <)2 dt 


Vi 2(At)^ 

Uttt { t , x ) 


'tk 


2 (At)" 


'4+1 


2At 


{tk+i - t) dt + ^ - t) dt. 


^fc+1 _ -fc 


From (lOl) . we obtain 
(6.5) 

Multiplying (16.511 by /r(4>''+^ — 4)'') and integrating over [0,L] yields: 


g'- -(- (j)fc 

AtT{’ = 


fL gfe+1 _ gfe 

h- 


At 


($fc+i _ dx 


\ [ /i ($'=+!)" dx- i dx- At [ ($'=+1 - $'=) dx. 

2 Jo 2 Jq Jo 


( 6 . 6 ) 


Furthermore, from (13.31) with t = tj,, 1 and Taylor’s Theorem, we get \/w G Hq{0,L): 


Ut(tk+l,x) -Ut{tk,x) 

u--- w dx 

^ At 


A '^xxifk+liX^ -\~ Uxxifk^X^ 

A --- dx 


L) Ut{ti;,L) Utxitk+l^L') UixitkiL) 

M -- w[L) + J -- Wx{L) 


At 


j Ux{tk+1, L) -\-Ux{tk, L) uifk^i, L) u{tk, L) 

+ ki - - - Wx { L ) + K 2 -^- W [ L ) 

I '^tx (tfc +17 L) “t“ '^tx ( 47 L) / X \ 1 (4‘+i; L) “t“ (4 7 L) / 7 - \ 

+ di- - - Wx{L) + d2 - - - w{L) 

Cl(4+l) + Cl(^fe) I C2(4 +i) + C2(4) ,x\ \4.rrkr \ 

+ Cl ---- Wx [L) + C2 -r- w{L) = At T 2 [w ), 


(6.7) 



























BEAM STABILITY AND SIMULATION 


25 


with the functional : Hq{ 0 , L) —>■ R defined as 
= 


rL I r*'=+i 


2{At) 


2 {tk+i-tf dt +dt\wdx 


A I _ t) _ M 


-\- k2 


2At 


+ M 


+ J 



'tk 


2At 


2(At) 


2(At)2 


2 (4+1 -tf dt + I 


(4+1 - dt+ [ 
Jt, 


2(At)2 

*“+2 Uttttx{t, L) 2 Jj . I / N | 

(4- - t) di I Wx(i) 


2At 


(4+1 -t) dt- f 
Jt, 


*''+5 Uttx{t,L) 




2At 


2A< 

(4+1 -t) dt- -t) dt \ w{L) 


(4 - t) dt Wa:(d2) 


+ di 


Utttxit,L) 

, , 2At 

''=+3 


(tfc+i -t) dt - [ 

J ti 


''+5 Utttxit,L) 

2At 

''+5 Uttt{t,L) 


(4 - t) dt i(;a;(L) 


r^fc + l 


+ Cl • 


2At 


(Ci)«(0 


+ C 2 • 


. , 2At 

■fc+i 

(C2)»(0 

f , 2At 


(tfc - t) dt i(;(L) 


<‘‘+' - - /( 2At 

(tfc+1 -t) dt- ^ ~ 


(tfc+i - t) dt - / ^ - t) dt I u;(L). 




2At 


( 6 . 8 ) 


Now, from (13.221) and (16.71) follows Vw/j G Wh- 


(6.9) 


Jo — 1^hdx + /g A *""" 2^""" {Wh)xxdx 

+M ^^"("i7^'°W (n;.)(£) + j£i^iW^(u;,),(L) 

+J:+i(L)+$^(L) , 1 /fN , . <I>*’+i(L)+4>*’(L) /r\ 

+di " ^ —^^(w/i)a;(L) + d2-- ^—^WhiL) 

+ci • + C2 • P^’^ WhiL) 


= -AtT^{wh) + G^Kwh), 

where the functional G\{wh) is given by 

G\{wn) := ^x 

(wh)x(L) 


( 6 . 10 ) 
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A Taylor expansion of Q about yields with (13.511 : 


( 6 . 11 ) 


Cl(*fc + l)~Cl(*fc) A Cl(^fc + l) + Cl(^A:) t '^tx {ik + 1 ,L)+Utx{tk iL) A + 'T'k 

At 2 2 ^^-^3 ’ 

C2(tfc + iy—C2(^fc) C2{tk + l)-\-C2{tk) Ut{tk + l,L)^Ut{tk,L) /\tT^ 


At 


with 


= 


r*fe+i /■‘'=+i (Ci)ttt(O 


Tf = 



(tfc+i - dt+ f 
J ti 


t, 2(At) 

*''+i {Ci)tt{t) 


(tfc - dt 


(tfc+i -t) dt- " %At ~ 

/■b 

( 4+1 -t) dt- j 

Jtk 


2 Utttx {t-i dj) 

2At 


{tk -t) dt\ , 


(tfc+i - 1 )^ dt+ f 

Jtk 


'=+1 (C2)ttt(0 ^i2 

dt 


(4+1 ^ ^2At ***' ^ 

fh 

(tk +1 -t) dt- i 

Jtk 


*'=+i Utility L) 


{tk -t) dt\ . 


Using (I3.23L (13.241) . and (16.111) . we get 


( 6 . 12 ) 


1, C.t^+cJ,! u 3>^+Ui)+‘f>4i) 

-Ai-- 2 -- 2 - 

1, C,t"+Cj,2 , 3.'=+UL)+3.'=(A) 

-At- ^2 - 2 - °2 - 2 - 


-AtTi-Gl 

-AtTf - G^, 


with 


G2" 

G3' 


4 utx{tk+itL) + uf^{tk,L) 
Ol -T- 


^2 


uf{tk+i,L) + u^{tk,L) 


2 
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In (16.911 we now take Wh ■= At 




G Wh, due to (16.5L Using (16.6|) and (j6.12p . yields: 


, fe+l||2_|| fc,|2 


1 At 

||2 = -(Atf 2 y A(e^f+ + + 


- (At)2 fci 


/o 

1 


(A) + es(^) 


-fc+l 


(Tf),(L) + fc2- 


(L)+e'=(L) . 


Tf(U) 


At / Ce,T^+Ce^ , . 

- - 2 -- 2 - j 

_ +^1 +Ce^ 

^ V 2 ) 2 2 ^2 


- Pi 


/•fe+1 I /"k 
Se,l "T Se,l 




At / 

T 


Ce.t^+Ce.2 , J $'=+l(L)+$'=(L)\ 

---h O 2 - 


- At<5s 


/$fe+i(L) + $fc(L) 




2 / 

Ce.t' + Ce" 


_ ^^£2 '^e,2 ~l~ ‘^6,2 p Ce,2^ + Ce,2 


- P 


^fe +1 J_ /-fc 

Se.2 S, 


e,2 




Therefore, 




I fe+l||2_|| fc,|2 


1 Af 

f < -(Atf-y^ A(4+i + eL)(Tf)..dx + —Gj(<i>'=+i + il>'=) 


_ (At)2 ^ 

aA;+ 1 _j_ 

- p^^A_^.((At)2T3'= + AtG^-) 




- P 


/■fc+1 I /•k 
Se,2 ^e,2 


• {{AtfT^ + AtG'^) 


(6.13) 

Next, from (I6.10p follows: 


- -(At)2T2''($''+i+ $'=). 


|G5;'($'=+i+$^')| < G ^ p ^fc||2^ 

|Ut(tfc+i,L) - M?(tfc,P) 2 , |UL(P+u A) - I 


At 


At 


{tfc+i,L)+u^,j,(tfe,L) 2 , |Ut(tfc+i,P)+ u^(tfe,L) 


-r + 


(6.14) 


' 2 ' ' 2 
+ |$'=+i(L) + $'=(L)|2 + |4>^+i(L) + $^(L)P 


< C (11$'=+! + + |$'=+i(L) + $'=(L)|2 + |$^+i(L) + $^(T)|^ 

I pifc + 1 

+ ■^ / Il'“tt(^)lli2 + |wtt(^) A)P + dt + ). 


'tk 


(6.15) 
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It can easily be seen that 


(6.16) 


|Tf < At 


C^k+1 


'ik 


pik + l 


\\uttt{t)\\H 2 dt < CAt / \\utttit)\\H 2 dt, 


'ik 


ftk + 1 

(6.17) \\T,Y<CAt / ||«,,,(t)||^, + ||(Ci)«|P + ||(Ci)m|pdt, 

Jtu 


(6.18) 

and 


'tk 
ptk + l 


P^k + 1 

WT^r < cAt / wuutmj,^ + ii(C2)«f + imutw^dt, 

Jtk 


r2'=($'=) < C\^\\^^\\i2 + \<s?^w + miw + 

rtk+i 

+ At / I|wtt(t)||ff4 + ||uttt(t)||^2 + ||Mtttt(t)||ff2 C?t 

Jtk 

/ tk+1 \ 

||(Cl)«(0f +ll(C2)«(t)lprftj. 

For the above estimate, we rewrote the second term of in (16.81) as: 


(6.19) 


Uttxx{t,x) 


2At 


/, ,, /* 2 Uttxxit^x) , , I ^ 

{tk +1 -t) dt- - — - {tk -t) dtj dx 


fifc+l 




2At 


'^ttxxxi^^ (Z/) + J' Uttxxxx{ij dx^ dt 

i'^ttxx{t.,L)^^(^L/^ Uttxxx{t.; J Uttxxxxit, x')^ dx \ dt^ 

using $^(0) = d)J;(0) = 0, and then the Sobolev embedding Theorem. From ()6.13l) - ()6.19ll . 
now follows: 


< C At{\\z':+^^ + \\z^J^) + At\\ut\\l 


ik + 1 


t\\C([tk,tk + lhH^) 


\\Uttit)\\L^ + WuitX)\ +\Uttx(tX)\ dt 


ik 


,2 ptk + i 

+ / ll(Cj)ttlP + IKCOtttlP 

i=l dtk 


pik + l 


{Atr + \\uut{t)rH^ + wntumr^ dt . 


'tk 


( 6 . 20 ) 

Let now nG{l,...,5'}. Assuming At < ^ (with C from (16.201) '). and summing (16.201) over 
e {0,..., n}, gives: 


+ C f At^ ll^ell^ + ll“t llZ[0.T];H2)+ ||Utt |||2(o 


2 


k=l 


+ (At)"^ ll(Ci)tt(0llz.2(0 y.gnl + ||(Ci)ttt(0llL2(0 -r.RnI 

J=1 

+ \Wtt{t)\\'i2(^Qx-m) + lkiti(^)lli2(0.T;Lr2) + l|wtMt(t)||i2(o,T;Lr2) 


( 6 . 21 ) 
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Finally, using the discrete-in-time Gronwall inequality and (EH), we obtain: 

||^n+l||2 ^ C* MIZgiP + 


+ ll(Cl)tt(OllL2(o r-Rni + ||(Ci)tti(OllL2(0.r:R^~) 

.i=l 

+ l|wti(Olli2(o,T;ffi) + ll^t«(OllL2(0,T;Lr2) + llL2(0,T;Lr2) 


( 6 . 22 ) 

The result now follows from (I6.22L (13.2511 . and the triangle inequality. 


□ 
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